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dS • Abstract 

o 

■ We have analyzed the effect of electron-LO phonon interaction in a double- 

ON ■ barrier resonant tunneling structure under a magnetic field B applied paral- 



lel to the tunneling current. While the low temperature anti-crossing phe- 



C3 

a 

*^ • nomenon has already been investigated, here we study the phonon absorption 

resonant magnetotunneling at finite temperatures. The Matsubara technique 



is used to sum up resonant diagrams of higher orders, and the result is selfcon- 
sistently renormalized. The phonon absorption tunneling spectrum has been 
calculated numerically. The result shows that the phonon-absorption process 
produces two broad inelastic wings on the main elastic tunneling peak, the 
strength of which is sensitive to the resonant condition, and grows with in- 
creasing temperature. The width of the inelastic wings are proportional to 
B 1 / 2 . This is in contrast to the appearance of anti-crossing in the phonon 
emission tunneling process. The difference is due to the fact that the phonon 
emission is a coherent process, while the absorption is a typical incoherent 
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process since absorbed real phonons have random phases. 
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I. INTRODUCTION 




re- 



sistance by Sollner et a/.EI, there has been increasing interest in resonant tunneling through a 
double-barrier resonant tunneling structure (DBRTS) both experimentally and theoretically. 
Among the important tools of characterizing DBRTS, one is the DC transporti~@ together 



the DC current noise of a DBRTS has also been investigated. A common feature of various 
resonant tunneling structures is the inelastic resonant tunneling due to the coupling of tun- 
neling electrons to phonons. This phenomenon is not only interesting but also useful because 
it provides a way to investigate both electron modes and phonon modes. A theory of phonon 
assisted resonant tunneling through a DBRTS in the absence of an external magnetic field 
was presented by Wingreen et a/.0 and by Jonson0. Assuming a perfect sample without 
imperfection and impurity, in these works an exact solution was obtained when the system is 
reduced to one-dimension and when tunneling processes through different quasibound states 
in the well are decoupled. The problem becomes much complicated when these tunneling 
processes are coupled, for example by phonons. 

Magnetic field is another important tool for sample characterization because of the for- 
mation of Landau level spectrum and the drastic modification of electronic wave functions, 
especially under strong external field. For a DBRTS, the challenging situation is when the 
magnetic B is applied parallel to the tunneling current I, as schematically illustrated in 
Fig. 1. To clarify the important features of Fig. 1, let us consider the experimentally most 
extensively studied DBRTS GaAs + /Al 03 Gao/iAs/GaAs/Al 03 GaojAs/GaAs + , with both 
the barrier width and the well width of the order of 40-60 A. In this case the barrier height 
is about 300 meV, and in the well there exists only one quasibound state with energy q,. 
We will set the zero reference energy at the bottom of the conduction band of the collector 
(labeled by c), and so the conduction band minimum of the emitter (labeled by e) is raised 
to eV by the bias V. Under a magnetic field B perpendicular to interfaces which are parallel 



with phonon induced or laser induced tunnelin; 



I. In connection to such characterization, 



3 



to the :n/-plane, the quasibound energy level e& splits into a series of discrete Landau levels 
€b + (n + ^)Huj c , where ou c = eB/m*c is the cyclotron frequency in GaAs. The electronic 
energy levels in the emitter and collector, on the other hand, form sets of Landau bands 
eV + e z (k ejZ ) + (n e + \)Tiuj c and e z (k C)Z ) + (n c + \)Tiuj c , respectively, where e z (k z ) = U 2 k 2 z /2m* 
is the kinetic energy along the tunneling direction. 

Such level splitting has important influences on both elastic and phonon assisted 
tunnelingllHl!. In a perfect DBRTS, if the electron-phonon interaction is ignored the tun- 
neling tunneling process conserves the Landau level quantum number (n e = n = n c ). If the 
magnetic field is sufficiently strong such that hu c is much larger than the thermal energy 
k-^T, the resonant tunneling is indeed one- dimensional with all tunneling channels decoupled 
from each other. The electron-phonon interaction mixes the channels. If the magnetic field 
is tuned into the resonant condition hu c = huo, where u Q is the optical phonon frequency 
in GaAs, then the resonant tunneling processes through two adjacent Landau levels are 
resonantly coupled. In fact, this is one way to observe the strong renormalization of Landau 
levels due to the electron-phonon interaction0'0. In a 3D electron gas, the renormalization 
effect is proportional to g 4//3 , while in a 2D system it is proportional to g, where g is the 
electron-phonon coupling constant to be defined later. 

Under a strong magnetic field and at low temperatures, one relevant interesting phe- 
nomenon is the formation of magnetopolaron. Let us consider the simpler situation that 
only the lowest Landau band in the emitter (n e = 0) is occupied with quasi Fermi energy 
e ej F- In a perfect DBRTS, one might expect to see in the I-V characteristics the elastic 
resonant tunneling peak around the bias eV — — e e F, with series of phonon replicas at 
eV = tb + nhu c — e e) F + vTiujq, where v is an integer. However, complication arises when 
the condition u c = iy — u')u; /(n' — n) is satisfied. Then, an electron in the emitter can 
tunnel resonantly into two coherent degenerate states in the well: into the Landau level n 
by emitting v phonons, or into the Laudau level n! by emitting v 1 phonons. The two final 
states are strongly hybridized because of the coherence, and so in phonon emission resonant 
magnetotunneling experiments, magnetopolarons have been observed^S A standard way 
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to present such experimental data is to plot resonant tunnelings peak position as a function 
of the magnetic field and the bias voltage. In this two-dimensional plot, the existence of 
magnetopolarons is reflected by the anti-crossing. 

A phonon emission resonant magnetotunneling process leading to two coherent final 
states, indicated as EO and El, is illustrated in part (a) of Fig. 2, where the thick lines form 
the Landau fan in the emitter, and the thin lines form the Landau fan in the well. Under 
this resonant condition, the I-V characteristics will exhibit two resolved peaks of almost 
equal strengthBi!. A similar diagram for phonon absorption resonant magnetotunneling is 
demonstrated in part (b) of Fig. 2. Here again we assume that only the lowest Landau level 
in the emitter is occupied, and the bias is tuned into elastic resonant tunneling eV = e^ — £ 6j f- 
Then, the two sets of Landau fans coincide with each other. At the magnetic field strength 
hu c = HoOq, an electron in the emitter having the proper energy can tunnel resonantly into 
the state AO elastically, or into the state Al by absorbing a phonon. We are interested in 
finding out how the energy levels AO and Al would be modified when they are resonantly 
coupled by phonons. In other words, we would like to derive the resulting transmission 
spectrum when the phonon absorption resonant magnetotunneling takes place. 

The essential issue here is the coherence. Since absorbed thermal phonons have random 
phases, they can not participate in coherent process. This problem was investigated in our 
recent work!!, where relevant references are listed. However, this work is based on a low order 
perturbation expansion, together with an exponential renormalization (see e. g. Ref. p3| ). 
In the vicinity of resonance, low order perturbation expansion is not sufficiently accurate. 
In the present work we improve our previous results with an infinite summation of resonant 
diagrams and a selfconsistent renormalization. Analytical expressions can be derived under 
the reasonable assumption of a dispersionless longitudinal optical (LO) phonon branch. 
To present a complete analysis, we will treat both the phonon emission and the phonon 
absorption resonant magnetotunneling. 

Sec. II describes our model Hamiltonian, from which the transmission probability and the 
differential conductance will be derived in Sec. III. The self-energy and the vertex correction 
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will first be analyzed at low order for the general case, and then be summed up to all 
orders for a model two-level system, which contains the essential physics for phonon assisted 
resonant magnetotunneling through a DBRTS. This will be done in Sec. VI, and the so- 
obtained self-energy is renormalized selfconsistently in Sec. V. Numerical results will be 
presented in Sec. VI, followed by some remarks in the last Sec. VII. The electron-phonon 
scattering matrix elements are calculated in Appendix A. 



To describe our model system, let us refer to Fig. 1 which specifies the potential and 
various energy levels of a DBRTS under a bias V, with zero reference energy set at the 
conduction band minimum of the collector. We assume the same material for the emitter, 
the well, and the collector, and also assume only one quasibound level et, in the well, which 
is usually the case for realistic samples. In our model we will ignore the scattering due to 
the interface roughness and impurities. Since the magnetic field B is applied in the positive 
direction of z-axis perpendicular to interfaces, it is convenient to choose the Landau gauge 
with the vector potential A=(0, Bx, 0). 

The eigenstates in the well are labeled by the Landau level quantum number n and 
the index of degeneracy k y . However, in the emitter (or collector), besides the similar 
quantum numbers (n e , k e ^ y ) [or (n c , k C;y )], we need also to specify the kinetic energy e z (k e>z ) 
[or e z (k CtZ )] of an electron along the tunneling direction. For the convenience of presenting 
mathematical formulas, we define the simplified notation a = (n, k y ) for the states in the 
well, and (5 = (rij, kj tV , kj tZ ), where j=e for the emitter, and j=c for the collector. In terms 
of these notations, the normalized eigenstates and the eigenvalues are 



II. THE MODEL 




(1) 




(2) 



in the well, and 
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(f)j,p(r) = exp(ikj, z z + ikj, y y)<p n (x + l 2 k jty ) , (3) 
= eV5 j>e + e z (k jtZ ) + fuo c (nj + -) (4) 



in the emitter (j = e) under the bias eV, and the collector (j = c). Here / = ych/eB is the 
magnetic length. 

Using this set of basis functions, our model system in Fig. 1 can be described with the 
tunneling Hamiltonian 

H = H c + H p h + He-ph . (5) 
The electronic part of the Hamiltonian with an external magnetic field is given by 

+ [^,^4cj,/9 + h.c] . (6) 

The tunnel matrix elements Vj } p a have to be calculated using the eigenstates listed above. 
Since the interfaces are assumed to be perfect, the momentum ky parallel to interfaces 
remains invariant during the tunneling through a potential barrier. Both the Landau level 
index and k y are then conserved, and so the calculation of the matrix elements Vj^ a reduces 
to the solving of a one-dimensional Schrodinger equation^, following with the application 
of the Bardeen's prescriptionil. 

Electrons interact with LO phonons. Here we will neglect the weak phonon dispersion, 
and the resulting phonon Hamiltonian is simply 

ftp h = hu frqfrq • ( 7 ) 

q 

Since electrons in the well has a finite life time, it is reasonable to expect that the dominating 
contribution to electron-LO phonon interaction comes from the situation when electrons 
occupy the quasibound state in the we 110. In terms of the Frolich Hamiltonianil, the 
electron-LO phonon interaction can be expressed as 



H c . ph = £ M QQ1 (q)(6t +6 q ) C t iCa , (8) 

where a = (n, k y ), a\ = (ni, /c y — q y ), and 

M Qai (q) = JL-J ^(rje^r)*, 

M 2 = 7re 2 ^ f — - -) . (9) 

The above Hamiltonian describes the general case. We are interested in the phonon 
absorption resonant magnetotunneling process as illustrated by part (b) in Fig. 2. Since 
the cross section for multiphonon absorption is much smaller than that for single phonon 
absorption, our problem can be well represented by retaining only the two lowest Landau 
levels n = and n = 1. With low impurity concentration in the emitter such that only 
the lowest Landau level n e = in the emitter is occupied, the resulting two-level system 
can be solved exactly. This exact solution can be extended to study the phonon emission 
resonant magnetotunneling if the process demonstrated by part (a) in Fig. 2 will be modified 
as follows. We first increase the impurity concentration in the emitter such that under very 
strong magnetic field, only the two lowest Landau levels n e = and n e = 1 are occupied. 
Next we set the bias at eV = e b + hu — e e) p. Then, when the magnetic field is tuned to 
Tiuj c = hu , electrons in the second Landau level (n e = 1) in the emitter can tunnel resonantly 
either into the El state via an elastic process, and into the EO state by emitting a phonon. 
This is again a resonant tunneling in a two-level system with two coherent degenerate final 
states. 

In the rest of this paper, we will first perform a general theoretical analysis. Then we 
will continue to study in details this exactly solvable two-level system. 

III. TRANSMISSION PROBABILITY 

To calculate the conductance of a DBRTS, we must first calculate the transmission 
probability of an electron from the initial state (n e , k etV , k £tZ ) in the emitter with energy e to 
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the final state (n c ,k c , y ,k CtZ ) in the collector with energy E\. This transmission probability 
can be expressed in terms of the Fourier transform of the two-particle Green's function as0 

T aai (e,ei) = 7 e (e)7 c (ei) 

x J d JI^S^^-^K aai (r, s, t) , (10) 

where 

K aai (r,s,t) = e(s)G(t) 

x(c Q (r- S )i(r)c ai (t)4(0)) (11) 

is the so-called transmission Green's function. Here 

(ni, = (n c , fc C]3/ ) label two electronic states in the well. 7_/(e) = 2^X^/3 | V} i( 3 Q | 2 <5(£ — -Ej,/?) 
is the escape rate of an electron from the well to the emitter j = e or the collector j = c. 
The total escape rate from the well is simply 7(e) = 7 e (e) + J c ( £ )- 

It can been shown0'0il that the Fourier transform K aai (e, Ei) of K aon (r, s,t) is in fact 
equivalent to a two-particle Green's function with input energies (e, E\ — Q) and output 
energies {e\, e — ft), analytically continued into the upper half plane of the complex variables 
(e, £1) and into the lower half plane of the complex variables (e — Q, E\ — Q). This two- 
particle Green's function is schematically illustrated in Fig. 3. The diagrammatic technique 
for an interacting electron-phonon system under a uniform external magnetic field has been 
extensively discussed in the literature§HH§, and here we will follow the conventional approach. 
After completing such calculation and then setting Q —>■ 0, we obtain K aai (e,ei) and 

T aai (e, £1) = le(£)lc(ei)K aai (e, £1) . (12) 

If there is no electron-LO phonon scattering, K aai (e,Ei) can be readily derived from the 
electronic Hamiltonian T~C e as 

iW^j- {£ _ Ea)2 + {l/2)2 - W 
The total transparency is defined as 



(14) 



with the corresponding correlation function 



(15) 



In the wide band approximation, where the energy dependence of 7 e and 7 C can be neglected, 
T a (e) and K a (e) are simply connected as 



Therefore, our task is to investigate the correlation function K a (e). 

In terms of diagrams, K a (e) is represented in Fig. 4. In this figure the arrowed double line 
towards right (or left) is the dressed one-electron retarded (or advanced) Green's function 
G\r(«, e) [or Ga(cx, e — Q) |n->-o] of the entire system including the electron-phonon interaction, 
where the horizontal arrowed single lines are the corresponding bare one-electron Green's 
functions G^(a,e) and G^\a,e — fi)| n ->o of the electronic system without the electron- 
phonon interaction. If an arrowed single line is tilted away from the horizontal direction, 
the variables in the corresponding Green's function should be summed up. All retarded (or 
advanced) Green's function are analytically continued into the upper (or lower) half plane 
of the complex variable e. Based on this diagram, K a {e) consists of two terms 



The first term is just a product of the two analytically continued dressed one-electron Green's 
functions 



T a (e) = 7e7c-Ka(e) • 



(16) 




(17) 




(18) 



The second term is the vertex part, which can be expressed as 




(19) 
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Let us first consider the Green's function Gs(a,e). To study the effect of electron-LO 
phonon interaction on the Green's function, we will use the Matsubara techniqueHS From 
the Dyson equation for the Matsubara Green's function 

£ (0) («,£) 



G{oc,e) 



where 5 = 7sign(Im[e]), and 

with discrete Fermion frequencies e = iir(2n + 1)/cbT, it is clear that we need to calculate 
the self-energy £ (a,e). While the self-energy will be analyzed in details later to all orders 
of electron-phonon interaction, to the lowest order shown in Fig. 5a, we have 

Z(a,E) = -k B T J2 g (0 Ha u s 1 )V ( - \c l ,e-s 1 ) 

ei,ai,q 

x|M Qai (q)| 2 , (22) 

where 

P"M= (23) 

with discrete Boson frequencies u = mlnk^T . 

To the same order as the above self-energy, the vertex correction JC(a,e) in (fK3|), with 
the corresponding diagram Fig. 5b, can be expressed as 



K{a,e) = 2nk B T £ | M aai (q) 



|2 



ot.\,e\,q 

xP (0) (£ - s 1 )g ( -°\a 1 ,e 1 )g^(a 1 ,e 1 - Q) . (24) 

The factor l-nk^T in front of this expression is due to the fact that the E\ integration in ( |Hj) 
has been replaced by a summation over the discrete Matsubara energies. Comparing this 
expression to (^2|) for self-energy, we expect a relation between JC(a,e) and S(a,e). This 
will be clarified in the next section. 
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Now we have set up the scheme to calculate the transmission probability with (f!6|)-(|T9|) 
and (pip. Knowing the transmission probability, the tunneling current through a DBRTS 
can be calculated from0 

x[f e (e-eV)-f e (e)]de, (25) 

where 

e - e j>F 



is the Fermi distribution function in the emitter (or collector) for j = e (or j = c) with the 
Fermi energy e e ^ (or e Ci p)- For realistic samples of DBRTS, resonant tunneling occurs at a 
bias eV much larger than the thermal energy k&T. Under this condition, we have f c {z) — 0. 
Hence, the differential conductance G(eV) = (dJ/dV) has the simple form 

«cv)-^jT ? w.,(-2&^Q)*. m 

IV. SELF-ENERGY AND VERTEX CORRECTION 



We will first calculate the lower order self-energy (|22|) and vertex correction (p4|) , and 
then generalize the results to the sum over terms of all higher orders. 

Let us start with the self-energy. The analytical continuation of the self-energy (22) has 
been performed by Abrikosov et.alS in details. The key step is to treat E\ as a continuous 
complex variable and rewrite the summation over e% into an integral over the path C as 
indicated by the solid arrowed curves in Fig. 6a. In this contour integral, the contribution 
comes from the horizontal path Im[ei] = and Im[ei] = e. The resulting self-energy is then 
analytically continued into the upper half plane of the complex variable e. We then obtain 

T, R (a,e) = 



27 



x Im 



-{e l^x(q) I 2 / 

Gf ef)] Df (q, e- ei ) tanh } 



+ Gf(a 1 ,e-e 1 )lm D^>(q,e 1 ) coth ( -f±- ) , (27) 



( £l 



\2k B T 
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where (q, e) is a bare phonon Green's function. The contribution from the pole at 
Si = is in Fig. 6a is included in the last term at the right hand side of (|27|) . Using the 
spectral representation 



1 

7T 



der 



Im 



Gft (a, Ei] 



Ei — £ — iS 



(28) 



and a similar form for /^(q, e), we arrive at the final form for the lower order self-energy 

E R (a,e) = — - V | M actl (q) | 2 [ du [ de x 
Im [D%\q,uj 



Im 



G R (a\, E\ i 



UO -\- E\ — E — i§ 
x [N(u) + 1 - /( £l )] , 



(29) 



where N(u) is the Planck distribution function. 

To further analyze the above expression, we will adopt a dispersionless bulk LO-phonon 
spectrum u;o(q) — wo, which is the commonly used approximation. In this case the phonon 
Green's functions becomes q-independent D^\q, u) c± D^\oj). Since the system is spatially 
homogeneous in the (x, y)-plane,the electronic Green's function does not depend on k\ y , so 
Gjt(ai,Ei) = Gfi(ni,ei). Consequently, in (p9|) the summation over (q, k ly ) acts on the 
matrix elements M aai (q) alone, and the result derived in Appendix |A] is 

M 2 



E I^Wq) \ 2 ~- 

fei a ,q 



(30) 



A nni = V27r^| nd(C 2(pi - p2) e- 
Pil Jo 



X 



(31) 



Here P2 {( 2 ) is the Laguerre polynomial, and p\ (or p-z) is the larger (or smaller) one of 
(n,ni). 

What left to be done in (p9| ) is the two integrations. The integration over uj is trivial 
because of the dispersionless phonon spectrum lm[D^ (u)} = —it [5{uo — uj ) — 5(u + uj )]. 
Since in realistic DBRTS samples the 7 is less than 0.5 meV which is much less than the 
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thermal energy k^T in the problem of our interest here, the integration over e\ can be carried 
out with the residual method. Together with the equality N(— uj) = — 1 — N(u), we obtain 
the simple form 



E fl (a,e) 



M 2 
4^1 



in 



N(hco ) + f(E ni ) 
e - E m + huj + 27/2 



+ 



iV(^o) + l-/(E rai ) 
e - E ni - hu + ry/2 

The important quantity is the imaginary part 

M 2 



(32) 



Im[E fl (a,e)] 



4tt 2 / 



7 
2 



N(hoj Q )+f(E ni ] 



(e - £ ni + ^ ) 2 + ( 7 2 /4) 
AT(^o) + 1 - 



(33) 



( £ -E ni -^ ) 2 + (7 2 /4). 
which will be discussed later. 

Now we come to the vertex correction K(a,e) in ([24]). The procedure of analytical con- 
tinuation is similar to that for the self-energy £^(a,£), but using the contour of integration 
in Fig. 6b. After all algebraic manipulations, we obtain 



dei 
2^ 



Ptanh 



£1 



2k B T 



+ Qcoth 



£1 



2fc B T 



(34) 



where 



7> ee Df(e - £l - O) Gi?(ni, £1 + ft) Im^m, £x)] 
+£>(«) ( e - £l - fi) ^(m^x) ImfG^K.e! + fi)] 

Q = Im[£>i?( £l )] 

xGg (ni, £ — £1) (m, £ — £1 — Q) . 

At this stage, one can set Q — > 0. Following the same argument which leads (^) to 
the vertex correction is derived as 



(35) 
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2 ™ ni 



+ 



e - E ni + ^ ) 2 + ( 7 2 /4) 
iV(^ ) + 1 - /( e - ^ ) 



(36) 



( £ _ jErii _^ )2 + ( 7 2 /4) _ 

Comparing this expression with (|33|), we see that at the lowest order of electron-phonon 
interaction, the self-energy and the vertex correction have similar structure. In particular, 
for the situation very close to resonance e ± u pa E ni , which is the main interest of the 
present work, we establish the following relation 

^, g ) = 47r|Im[ ^ K£)]l (37) 

7 

between the self-energy and the vertex correction. It is important to point out that the 
above relation (|37|) holds up to any order of electron-phonon coupling, provided that the 
phonon spectrum is assumed dispersionless0i3il. Consequently, from now on we only need 
to analyze the single electron self-energy. 

All above results were derived with the only approximation that the weak dispersion of 
the LO phonon spectrum can be ignored. If the tunneling is not resonant, these results 
are accurate enough to be compared with experiments at all temperatures. They are also 
reliable if at low temperature the thermal energy is much smaller than the LO phonon en- 
ergy. However, since we are interested in resonant magnetotunneling at finite temperatures, 
higher order electron-phonon interaction must be included in the self-energy and the vertex 
correction. This problem is exactly solvable if only two Landau levels n = and n = 1 in 
the well dominate the tunneling process. As described at the end of Sec. II for our model, 
this two-level system can be realized by fixing the concentration of impurities in the emitter. 
For the phonon absorption resonant magnetotunneling process, the tunneling electrons are 
emitted from the Landau level n e = in the emitter. On the other hand, for the phonon 
emission resonant magnetotunneling process, the emitting Landau level is n e = 1. From 
now on we will restrict ourselves to such two-level system. 

The procedure of calculating the contribution to self-energy by any order of electron- 
phonon interaction is the same as that by the lowest order. Analysis of analytical con- 
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tinuation leads to the following rule of constructing diagrams. We introduce two kinds of 
phonon lines: a dotted line describes a absorption process with a Bose factor N(huo), while 
a dashed one corresponds to a phonon emission with a factor N(huJo) + 1. Then, each dot- 
ted line should be followed by a dashed line, and vice versa. Consequently, all even order 
contributions vanish, and only the diagrams with odd number of phonon lines survive. As 
an example, the third order diagram is shown in Fig. 7. Based on the diagrammatic anal- 
ysis, the uth order contribution contains 2v electron-phonon vertices, v factors g 2 A§\ with 
g 2 (B) = M 2 /4ii 2 l, v phonon factors N(u Q ) + 1/2 =f 1/2, and 2u — 1 electron Green's func- 
tions. Therefore, the final form of each order's contribution is very similar to fl32|), except 
that there is only one term in the summation of n\, and Aq\ ~ 1. 

We will not present here the detailed calculation, but will only show the final results. 
Since our model neglects the electron-electron interaction in the well, it is equivalent to the 
case that in the well the occupation probability of each electronic eigenstate is negligible. 
Then, including all orders in g, the electron self-energy is obtained as 

T, R (a, e) = 

g2 (2s+l) N s+l( N+1 y 



E 

+ 



(e - E m + hcu + n72) 2s+1 (e — E n + 2 7 /2) 2s 
g 2(2s+i) N s( N + 1 y+i 



(38) 



(e - E ni - hcuo + * 7 /2) 2 * +1 (e-E n + z 7 /2) 2s 
where iV = N(huoo) is the Planck function. When e ~ Eq and so E ni = E±, the first term 
in (153) gives the phonon absorption resonant magnetotunneling. On the other hand, when 



e ~ E\ and so E ni = Eq, the phonon emission resonant magnetotunneling is represented by 
the second term in (|38T). 



V. SELF-CONSISTENT RENORMALIZATION 

Studies on resonance effects at very low temperatures in both free electron gas0@ and in 
DBRTSi-i indicate that resonant states are strongly coupled, and the perturbation expan- 
sion should be renormalized. A so-called self-consistent renormalization can be used, which 
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replaces all bare electron Green's functions in the self-energy formula with the corresponding 
dressed ones, resulting in self-energy equations. 

Let us first renormalize the first order self-energy to review the phonon emission process 
in which the tunneling electron starts from the Landau level n e = 1 in the emitter. The 
self-consistent equation is simply 

a 2 

SflK e) = r 7 n 7- ■ (39) 

Very near to the resonant e — E — Tiujq = and huo = hu c , this self-consistent equation has 
the solution 

E R (a,e)~(E + tku -e)/2±g, (40) 

which will produce the Landau levels anti-crossing phenomenonS'il. For a 2D free electron 
gas. the higher-order corrections has been analyzed in Ref. [32. 

Now let us turn to the main theme of the present work: finite temperature phonon 
assisted resonant magnetotunneling. Here both phonon emission and phonon absorption are 
involved, and we must use (|38f) as the basis for the renormalization. If we define £_r(ti, e) = 
£r(o!, e) and Gn(n,e) = Gr(qz,e) for a = (n,k y ), then the self-consistent equation can be 
expressed as 

oo 

e + tku ) = J2 9 2{2s+1) N s (hu; ) [N(hu ) + 1} S+1 G^ +1 (0, e)G%(l, e + hu ) , 

oo 

E fl (0, e) = £ g 2{2s+l) N s+1 (huj ) [N(hu ) + l] s G| s (0, e)G| s+1 (l, e + hco ) . (41) 



s=0 



Let us define the temperature dependent coupling constant // = g 2 J N {Tiujq)[N {Tiujq) + 1], 



the dimensionless thermal factor v = J N (huo) / [N (Tiujq) + 1], and 



le; -l- A i 2 [G Jl (O l e)G Jl (l,e + fiw b )] 2 " 1 J 



If we further define 

e- (E + E 1 -huj -ij)/2 



9(e) = 
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then, we can rewrite (fUD as 



E fl (l, e + hu ) = y/y[9{e) + 5} i + ^ e y u 

E fl (0, e ) = V?«[g(e)-g] 1 ^ ) (e) 
^ * \ 1 l + £(e)/f 



Substituting the last two equations of fl43|) into (|42l), we obtain the equation 



[1± Vl + 4L 2 ][l + L(z/+l/z/) + L 
2L 



21 



6 2 {e) - <5 2 (44) 



to solve for L(e). Knowing L(e), we go back to the first two equations of (|43f) to calculate 
the renormalized self-energies + frcjo) and E#(0, e). 



VI. RESULTS 

We will solve (|4"4] ) for a temperature corresponding to 7 C y^/J. There are several 
branches of solutions, and the proper one should be chosen according to the correct analytical 
properties and asymptotic behavior of the Green's functions. In general, the asymptotics 
of Green's functions at large e is 1/e. Consequently, the correct branch of solutions should 
correspond to lim^oo L(e) = 0. 

Let us first check a limiting resonant tunneling case. At | |^(^)| 2 — S 2 |<C 1, the proper 
solution corresponds to |1 + L(e)/v\ <C 1. Therefore, we obtain L(e) ~ —u, which yield 
the solutions Gr(1,e) — > and Gr(Q,e) — > ^/Jl(l — v 2 )/[6(e) — 5}. We see that in the 
phonon absorption process, the original elastic resonant tunneling peak in the transparency 
spectrum survives, and its position is not affected by the electron-phonon scattering. This 
is in contrast to the phonon emission process, where the original elastic resonant tunneling 
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peak splits symmetrically into a doublet of almost equal strength. At sufficiently high 
temperature k-^T ^> hu} , v — > 1 and then both G_r(0, e) and Gr(0, e) approach zero. 

This solution of limiting case will be used as the initial input when solving the selfcon- 
sistent equation ( |4*4"D numerically. We will study the case that one electron tunnels into 
the well from the n e = Landau level in the emitter, under the exact resonant condition 
6 = 0. We also choose 7 = 0.1 lg. Two temperatures are used in our calculation: T=174 K 
(corresponding to ^=0,3) with Im[6 l (e)]=0.1, and T=74 K (corresponding to z/=0,006) with 
Im[0(e)]=O.23. For a given value of Re[0(e)], we solve (fyj) for both Re[L(e)] and Im[L(e)], 
and then plot them in Fig. 8. The dotted contour is for T=174 K with Re[#(e)] varies along 
the contour from -4 to -0.2, and then from 0.2 to 4. For the other temperature T=74 K, the 
solution is the crossed curve with Re[0(e)] varies from -5 to -0.4, and then from 0.4 to 5. 

Based on the solution in Fig. 8, the corresponding dimensionless transmission probability 
K (e)g 2 [defined by (^) and (|T6|)1 is calculated, and the results are shown in Fig. 9a as 
functions of reduced energy Re[6 l (e)] JJl/g. In the absence of electron-phonon scattering, 
the sharp peak of transmission probability is marked as the curve 1. When the electron- 
phonon interaction is turned on, two inelastic wings are formed with increasing strength 
as the temperature is raised from T=74 K (curve 2) to T=174 K (curve 3). In order to 
demonstrate these inelastic wings more clearly, we use (|i~7|) to decompose the transmission 
peak K (e) into the main peak Kq(e) and the inelastic wings Kg(e). The results are plotted 
in Fig. 9b where curve 1) shows the main peak (for both temperatures), and where curve 
2) and curve 3) show respectively the inelastic wings for T = 7AK and T = Y1AK. For 
both temperatures, within the accuracy of our plotting, the broadening of Kq(e) due to the 
virtual phonon processes can not be detected because it can not be distinguished from the 
pure elastic peak. On the other hand, the temperature dependence of the inelastic wings is 
quite dramatic. 

The phonon absorption process is entirely different from the phonon emission process, 
which we have also calculated here and found the separation between the two resolved sharp 
peaks being proportional to 2g. The physical origin of this difference is that the phonon 
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emission is a coherent process for which the electron-phonon system remains coherent, while 
the phonon absorption is typically an incoherent process since absorbed real phonons have 
random phases. 

VII. REMARKS 

As demonstrated in the previous section, the absorption of real phonons during the 
resonant tunneling yields only minor corrections in the elastic main peak of the transparency 
spectrum. On the other hand, the inelastic process indeed produces two broad wings in 
the spectrum, the strength of which is sensitive to the resonant condition and grows with 
increasing temperature. If the tunneling current and so the differential conductance G(eV) 
can be measured in this way, we can further extract informations of the dynamical tunneling 
process as follows. Since the tunneling is dominated by the transmission probability of only 
one channel, T (e), fl26|) can be inverted to the form 

'2fi roo roo 

T (e) = — — di dfi 

e z K-Q± JO J-oo 

fe-e BtF -fM \ sinh£ 
x cos [ ^ ej —GM • (45) 

Substituting the measured G(eV) into the integrand and performing a numerical integration, 
the transparency spectrum To(e) can then be measured. 

The self-energy ([33]) has a similar form as that obtained by Fertig et.al.0 for a resonant 
tunneling system interacting with impurities under an applied magnetic field. While in both 
cases the magnetic field broadens the width of the resonant tunneling peak because of the 
violation of selection rules for the parallel momentum, there are two important differences. 
In our case the broadening reaches its maximum value when e = E ni ± Tiujq. Furthermore, 
the self-energy due to the electron-phonon interaction is proportional to B 1 ^ 2 , but for the 
electron-impurity scattering it is proportional to B. 

The present work is partially supported by the Norwegian Research Council, Grant No. 
100267/410. 

20 



APPENDIX A: SUMMATION OF ELECTRON-LO PHONON MATRIX 

ELEMENTS 



We will calculate the (ki y , q) summation of the electron-LO phonon matrix elements 

M 2 r 

I M aai (q) \ 2 = ^2 -— / dx dy dz dx' dy' dz' 



xe-Vev'r a ^)Mrmr')<t>«Ar') ■ (Al) 

First, we transform the summation into integrations. Since the Green's functions are inde- 
pendent of k y , we can conveniently choose k y = 0. Therefore, we have 

£|AW q )l^/|$/^|AW q )P 



fcij,,q 



= / 7~ — 7~ / dxdx'dzdz' 
J (2ir)- i q l J 

x <(^ - l 2 q y )^ ni {x' - l 2 q y )(Pn(x?)(p n (x). (A2) 



The q integration can now be performed to give 

c% exp [^(x - x') + - z')\ 



I 



(M 2 9f + ql + ql 



2tt 



q y | v/(z - ^') 2 + x') 2 



(A3) 



where X is the McDonald function. 

For realistic DBRTS samples, the width d of the quantum well is much smaller than the 
magnetic length /. In this case the difference z — z'm the argument of the McDonald function 
can be ignored, and the integration over z and z' is equal to 1 because of the normalization 
condition. By introducing two new dimensionless variables £ = q y l, r] = (x — x')/l, and 
( = x// we arrive at the final expression 

M 2 

£ I M Qm (q) | 2 = -^i nni , (A4) 

where 
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A nni = J d(d£dr)K (\ rit |) 

x«,(c-^)^(c-^) 

x«(c + ^)ft.(c + ^). (A5) 

Substituting the harmonic oscillator wave functions 

/ i \ 1/2 

into the above integrand, where H n (() is the Hermite polynomial, we obtain the results 
given by 
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FIGURES 



FIG. 1. A schematic illustration of a double-barrier resonant tunneling structure (DBRTS). 

FIG. 2. Diagrams illustrate phonon-assisted resonant magnetotunneling as the magnetic field is 
tuned into the resonant condition where the Landau level separation fiuj c equals the optical phonon 
energy Tiloq: A) A phonon emission process leading to two coherent final states (EO and -El), and 
b) a phonon absorption process leading to two final states (AO and ^41). 

FIG. 3. Diagrammatic representation for the two-particle transmission Green's function 
K aai (e, £i). 

FIG. 4. Diagrammatic representation for the total transmission K a (e) given by Eqs. (|17|)- 
©■ 

FIG. 5. Lowest order diagrams for (a) the self-energy E(a, e), and (b) the vertex Green's 
function K.(a, e). 

FIG. 6. Contours of integration for (a) the self-energy S(a,e), and (b) the vertex Green's 
function JC(a, e). 

FIG. 7. Resonant third order diagram for the self-energy. A dotted line denotes a phonon 
absorption process, while a dashed line denotes an emission process. The electron lines are marked 
by their Landau level numbers. 

FIG. 8. Solution Re[L(u)] and hn[L(co)} of Eq. (gj) for 5 = and 7 = 0.11c/. Dotted contour 
is for temperature T=174 K, and crossed contour is for T=74 K. 
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FIG. 9. Phonon absorption resonant magnetotunneling spectrum corresponding to the solution 
in Fig. 8. A) The dimensionless transmission probability Ko(e)g 2 is plotted as function of the 
reduced energy Re [0(e)] y/Ji/g. Curve 1) shows the transmission probability in the absence of 
electron-phonon scattering, while curve 2) and curve 3) show respectively the results at T = 74K 
and at T = 174K. B) Transmission spectrum decomposed into main peak K$(e)g 2 (curve 1) and 
inelastic wings K%(e)g 2 at T = 74K (curve 2) and T = 174k (curve 3). 
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